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ABSTRACT 

We present a unified model of infrared (IR), optical, ultraviolet (UV), and X-ray light curves for the 1983 
outburst of GQ Muscae (Nova Muscae 1983) and estimate its white dwarf (WD) mass. Based on an optically 
thick wind model of nova outbursts, we model the optical and IR light curves with free-free emission, and 
the UV 1455 A and supersoft X-ray light curves with blackbody emission. The best fit model that reproduces 
simultaneously the IR, optical, UV 1455 A, and supersoft X-ray observations is a 0.7 ± 0.05 M© WD for an 
assumed chemical composition of the envelope, X = 0.35-0.55, Xcno = 0.2-0.35, and Z = 0.02, by mass 
weight. The mass lost by the wind is estimated to be AM w j n d ~2x 10~ 5 M Q . We provide a new determination 
of the reddening, E(B-V) = 0.55 ±0.05, and of the distance, ~ 5 kpc. Finally, we discuss the strong UV flash 
that took place on JD 2,445,499 (151 days after the outburst). 

Subject headings: novae, cataclysmic variables — stars: individual (GQ Muscae) — stars: mass loss — ultra- 
violet: stars — white dwarfs — X-rays: binaries 



1. INTRODUCTION 

Despite their overall similarity, the optical light curves 
of classical novae show a wide varie ty of timescales and 
shapes (e.g., iPayne-G aposchkiril 119571) . Various empirical 
time-scaling laws have been proposed in the attempt to rec- 
ognize com mon patterns and to u nify the nova light curves. 
For example, McLaughlin] d 19421) proposed a compression of 
the time scale obtained by normalizing the time to t m , the 
time for the brightness to decrease by m magnitudes (usu- 
ally m = 2 or 3). The results were however unsatisfactory 
because slow novae reached their late stages relatively ear- 
lier than faster novae. A diffe rent approach was adopted by 
IVorontsov- Velvaminov I ([19483 who, based on a large collec- 
tion of nova light curves, found that m(t ), the visual magnitude 
at time t, is conveniently represented, on average, by m(t) = 
m.Q + b\ log(f-fo) at early stages and m(t) = nt\ + £>2 log(? - fo) 
at later stages with b\ = 2.5 and Z?2 > b\. 

The underlying problem of both the above approaches is 
that, during the decline phase from maximum, the flux in the 
visual band, initially dominated by free-free emission, be- 
comes more and more affected by the increasing contribu- 
tion from the emission lines, which finally dominates. The 
presence of these two competing and heterogeneous emission 
mechanisms is most likely the main cause for the difficulty 
to find a suitable normalization parameter for the nova light 
curves. A more appropriate way to monitor the evolution of 
the visual continuum is to use the Stromgren y filter, which is 
designed to avoi d strong emission lines and, in particular , the 
[OIII] line (e.g. jKalerlll98llLockwood & Millis Ifl976h . In 
the case of VI 668 Cvg. iKalen d 19861) has indeed shown that 
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the continuum y flux declines much faster than the H/3+fO III] 
fluxes. Consequently, the decline rate is much faster in the 
y band than in the visual band. This was demonstrated by 
iKaler I d 1986b . who showed that the y magnitude of VI 668 
Cyg declines with a slope of b\ =3.5 at early stages and it 
then changes to £>2 = 6.5 at lat er stages, values that are s ensibl y 
larger than those reported bv lVorontsov-Velyaminovl (f l948). 

The observations of VI 500 Cyg by Gallagher & Nev 
(1 19761) suggest that the visual and infrared continua are well 
represented by free-free emission d uring the early decay 
phase. The same applies to GQ Mus: Krautter et al. ( 1984) 
reported that, during the early 40 days (diffuse enhanced 
and Orion-phases), the energy spectrum of 0.3 -10// was 
well represented by optically thin free-free emission (see also 
iDinersteinl 1 1 9861) . The flux in the optical and infrared ap- 
pears to decay with time as F\ oc t~ a in various speed cla sses 
of novae (e.g.. lEnnis etlDl977l:IWoodward et al. II 19971 for 
V1500 Cyg). 

An interesting attempt to unify the light curves of novae 
was proposed by Rosenbush ( 1999a b), who suggested that 
the main parameter determining the shape of the light curve is 
the radius of the ejected shell, log7? s h e ii = log(f)+log(v), where 
t is the time from optical maximum, and v is the velocity of the 
ejecta. He plotted various nova light curves in the m-log/? S h e ii 
diagram and found that light curves overlapped each other. 
The rather large scatter of the data in that diagram is probably 
due to th e fact that the expansio n velocity of novae varies with 
ti me ( see Cassatella et al. 1120051) . 

ICassatella et al. I d2002l) studied the time evolution of the 
UV continuum flux in twelve CO and ONeMg novae and 
found a strong correlation between the tj, time and the time 
of maximum flux in the 1455 A continuum (see Fig. 4 i n their 
paper). In a subsequent paper, Cassatella et al. (2005) stud- 
ied the time evolution of the UV emission lines in seven CO 
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novae and found a strong correlation between the line ioniza- 
tion potential and the time of maximum emission normalized 
to the ?3 (Fig. 5 in their paper). These results strongly suggest 
that novae do indeed evolve following a common pattern be- 
ing time-normalized by the tj, time, whi ch mainly depends on 
the white dwarf mass (e.g..lLiyiolll992l) . 

Recently, Hac hisu & Kato I (|2006, hereafter referred as Pa- 
per I) found that the visual and IR light curves of several novae 
follow a universal decline law, and interpreted that in terms 
of free-free emission (Paper I; lHachisu & Kato I [20071 here- 
after referred as Paper II); in particular, the time-normalized 
light curves were found to be independent of the white dwarf 
mass, the chemical composition of ejecta, and wavelength. 
They also showed that the UV 1455 A light curve, interpreted 
as blackbody continuum, can also be time-normalized by the 
same factor as in the optical and IR. The authors, in the above 
quoted papers, determined the white dwarf mass and other 
parameters for a number of relatively well-observed classical 
novae. 

In the present paper, we apply the above universal decline 
law to the classical nova GQ Muscae 1983. GQ Mus is 
the first classical nova in which the superso ft X-ray turnoff 
was d e tected with the X-ray satellite ROSAT dQgelman et al. I 
fl99l lBalmanetal.1 [1991 lOrioetal. I I200i1k This ob- 
ject was also well o bserved by the UV satellite IUE (e.g., 
Krautter et all 119841) and in near-IR bands of JHKL (e.g., 
Whit elock etal. 111984 . 

The next section summarizes the basic observational char- 
acteristics of GQ Mus. In §3, we revisit the UV light curves of 
GQ Mus obtained with the IUE satellite. In §4, we briefly in- 
troduce our method based on our optically thick wind model. 
Light curve fittings of GQ Mus in the X-ray, UV, optical, and 
near IR bands are shown in §5. Discussion and conclusions 
follow in §§6 and 7, respectively. Finally, in the Appendix, we 
assess the problem of the contribution of the emission lines to 
the visual and infrared photometric bands for GQ Mus and for 
other well known novae: V1500 Cyg, V1668 Cyg, and V1974 

Cyg. 

2. BASIC OBSERVATIONAL SUMMARY OF GQ MUSCAE 1983 

GQ Mus was discovered by Lill er on 1983 January 18.1 4 
UT (JD 2 , 445,3 52.64) at my « 7.2 dLiller & Overbeekll983l) . 
iBateson I £l983) reported two pre-discovery magnitudes, i.e., 
my k, 7.6 obtained by Gainsford on 1983 January 15.597 UT 
(2,445,350.097) and my « 7.9 obta ined by Pattie on 1983 Jan- 
uary 15.628 UT (2.445.350.128). IKrautter eTaTl (fl98l sug- 
gested that the outburst took place 3-4- days before the dis- 
covery. In absence of precise estimates, we assume that the 
outburst took place at f B = JD 2,445,348.0 (1983 January 
13.5 UT), i.e., 4.6 days before the discovery by Liller on Jan- 
uary 18.14, and adopt f B = JD 2,445,348.0 as day zero in the 
following analysis. 

Early optical UBVRI and near infrared JHKL observations 
of GQ Mus show a fast evolution with ~ 18 days and 
t 3 = 48 days (IWhitelock et al. II1984I) . IKrautter et al~l (119841) 
observed GQ Mus over a wide wavelength range from 0.12 - 
10 [i during the early 40 days (diffuse enhanced and Orion- 
phases). These authors reported an extraordinary large am- 
plitude brightening by Amy ~ 14 mag, with respect to the 
prenova magnitude of my > 2 1 . After about 4 months of slow 
decline, the brightness stabilized in the range my ~ 10-11 
for about 400 days. After that, the nova gradually declined to 
my ~ 14 for another 400 days (see figures below). 

In this paper we will use the optical and IR photometric 



data from lWhitelock et al. I d 19841) and lKrautter et al. I d!984l) . 
which cover about 300 days after discovery. No other system- 
atic photometric observations are available at later stages ex- 
cept for the visual magnitude measurements (Vfes) obtained 
with the Fine Error Sensor (FES) monitor on board IUE, 
which cover 500 days after discovery, and the visual photo- 
metric data collected by RASNZ (Royal Astronomical Soci- 
ety of New Zealand) and by AAVSO (American Association 
of Variable Star Observers), which cover 1000 days after dis- 
c overy. 

iDiaz & Steiner I JT989) determined a photometric orbital 
period of 0.0594 days fo r GQ Mus, the shortes t known for 
a classical nova. Later on. lDiaz & Steiner] ( fl9 94) revised the 
ephemeris of GQ Mus based on the 1989 and 1990 photomet- 
ric observations, i.e., 

HJD min = 2,447,843.4721+0.0593650 XE. (1) 

We have adop t ed this eph emeris for our binary model. 
IDiaz & Steiner I d!989l 1 19941) also suggested that GQ Mus is 
a polar system because of the presence of X-ray emission and 
high (He II A4686)/H/3 ratio. 

Based on the optical observation from 1984 March 
to 1988 March (400 - 1800 days after the outburst), 
IKrautter & Williams I d 19891) suggested that the observed pro- 
gressive increase in the ionization level of the nova shell 
was due to a very hot radiation source with a temperature 
of > 4 x 10 5 K, increasing wi th time while the bolo metric 
luminosity remained constant. Ogelm an et al. I |l987) stud- 
ied the soft X-ray EXOSAT light curve from optical maxi- 
mum to about 900 days after it, and suggested a model of 
a very hot white dwarf remnant with a maximum temperature 
of (2-4) x 10 5 K evolving at constant bolometric luminosity. 

ROSAT observations clearly show that the supersoft 
X-ray flux decayed about 10 yr af t er the outburst 
(Ogelman et al. lll993HShanlev et al. ill 9951; lOrio et al. 11200 lb 
iBalman & Krautterll2001l) . This information is very impor- 
tant because, compared with our model predictions, it allows 
one to determine the white dwarf mass quite accurately (see, 
e.g., Paper I). 

3. UV OBSERVATIONS 

One of the most puzzling features of GQ Mus is its long 
lasting brightness in the UV range, which made possible its 
regular monitoring by IUE over as much as 11.5 yr. Because 
of satellite constraints, the IUE observations could not start 
before JD 2,445,385.07, i.e., 32.4 days after discovery. The 
early IUE observatio ns of GQ Mus are discussed in detail by 
Krautter et al. ( 1984). In the following we revisit the prob- 
lem of the color excess E(B-V) of GQ Mus and describe 
the long term evolution of the ultraviolet continuum and of 
the emission lines, as well as the secondary outburst which 
took place about 151 days after the main one. The ultraviolet 
spectra were retrieved from the IUE archive through the INES 
{IUE Newly Extracted Spectra) system 1 , which also provides 
details of the observations. The use of IUE INES data is par- 
ticularly important for the determination of the reddening cor- 
rection because of the implementation of upgraded spectral 
extraction and flux calibration procedures compared to previ- 
ously published UV spectra. 

3.1. The reddening correction 
1 http://ines.laeff.esa.es 
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GQ Muscae E(B-V) = 0.55 




2000 2500 
Wavelength (A) 

FIG. 1. — IUE spectra of GQ Mus obtained at four epochs, i.e., 108, 
151, 202, and 538 days after the outburst. The spectra have been corrected 
for reddening using E(B-V) = 0.55. The vertical dotted lines represent the 
wavelengths AA 1512, 1878 and 2386 A at which the extinction law takes 
the same value. With the adopted value of reddening, the stellar continuum 
underlying the many emission lines is well represented by a straight line all 
over the full spectral range except for the spectrum of day 538 above 2700 
A, due to the increased contribution from hydrogen Balmer continuum. Satu- 
rated points in the emission lines are denoted by pluses. Note the comparably 
harder spectrum at the time of the secondary outburst on day 151. 

The color excess of GQ M us has been determined from 
the hy drog en Balmer lines by Ide Freitas Pacheco & Codina I 
dl985l) and iPequignot et al. I (fl993l) . who found E(B- V) = 
0.43 and 0.50 ± 0.05, respective ly. About the same r ange 
of values has been reported by iRrautter et al. I (1 19841) and 
iHassallet al.ld 19901) . who found E(B-V) = 0.45 and 0.50, re- 
spectively, based on the strength of the dust absorption band 
around 2175 A in the early IUE spectra. Given the critical- 
ity of the reddening correction, we have attempted to improve 
its accuracy by using two independent methods based on the 
shape of the UV continuum and on the emission line intensi- 
ties. 

The Galactic extinction curve (ISeatonlll979h shows a pro- 
nounced broad maximum around 2175 A due to dust absorp- 
tion, but it takes the same value X(X) = A(X) / E(B - V) w 8 at 
A = 1512, 1878, and 2386 A, so that the slope of the straight 
line passing through the continuum points at these wave- 
lengths is insensitive to E(B-V) in a (A,logF(A)) plot. This 
circumstance can be used to get a reliable estimate of E(B-V) 
as that in which the stellar continuum becomes closely linear 
in the 15 12-2386 A region, and passes through the continuum 
points at the above wavelengths. From 8 pairs of short and 
long wavelength IUE spectra taken during the nebular phase 
we have in this way found E(B - V) = 0.58 ± 0.04. 

A different way to estimate the E(B - V) color excess 
consists in using the observed relative intensities of the 
He II 1640 A Balmer line, and the 2734 and 3203 A 
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FIG. 2. — Time evolution of the continuum fluxes at 1455 A and 2885 A, 
of the ultraviolet color index C(1455 — 2885), and of the visual flux. Fluxes 
are in units of 10~ 13 erg cm -2 s -1 A -I , not corrected for reddening. Only 
the color index has been corrected for reddening using E(B- V) = 0.55. The 
vertical dotted and dashed lines indicate the date of discovery and that of the 
secondary outburst, respectively. 

Paschen recom bination lines, and comp are these with the- 
oretical ratios (Hummer & St orevlll987l) . The intensity ra- 
tios /( 1 640) //(2734) and /( 1 640) //(3203) were measured in 8 
pairs of short and long wavelength IUE low resolution spectra 
obtained during the nebular phase. From a comparison with 
the theoretical values for an electron temperature and density 
of (T e , Ag=(20,000 K, 10 6 cm" 3 ) we have obtained 22 inde- 
pendent determinations of the color excess leading to a mean 
value of E(B-V) = 0.51 ±0.06. In the following we will adopt 
the error-weighted mean of the two above determinations, i.e. 
E(B-V) = 0.55 ±0.05. Examples of IUE spectra of GQ Mus 
corrected with E(B-V) = 0.55 are reported in Figure Q] for 
the following days (after the outburst): 108, 151 (secondary 
outburst; see later), 202 and 538. 

3.2. Evolution of the UV continuum 

We have measured the mean flux in two narrow bands 
20 A wide centered at 1455 A and 2855 A, selected to 
best represent the UV continuum because little affected 
by emission lines (Cassatella et al. 2002). Figure [2] 
shows the time evolution of the F(1455 A) and f\2885 
A) fluxes and of the UV color index C(1455-2885) = 
-2.5 log[F( 1455 A)/F(2885 A)]. The measurements were 
made on well exposed low resolution large aperture spectra. 
Unfortunately, the IUE spectrum of day 37 (SWP 19299), very 
important because it was the first obtained, was heavily satu- 
rated around 1455 A. To make a crude estimate of the corre- 
sponding flux we have determined the scaling factor between 
this and the next spectrum obtained on day 49 (SWP 19383), 
in a region where both were well exposed (1520-1620A), and 
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FIG. 3. — Time evolution of the observed fluxes in the most prominent 
permitted emission lines of GQ Mus, in units of 10 -12 erg cm -2 s . The 
vertical dotted and dashed lines indicate the date of discovery and that of the 
secondary outburst, respectively. 
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FIG. 4. — The figure shows the P Cygni profile of the Si IV doublet AA 
1393.74, 1402.77 A of GQ Mus at the time of the secondary outburst on 
day 150.64 (upper panel) and, for comparison, that of the Neon nova VI 974 
Cyg obtained on day 45 after discovery (lower panel): the IUE images used 
are SWP 20211 and SWP44390, respectively. The observed profiles have 
been fitt ed with theor etical P Cygni pro files obtained through the SEI method 
(Lamers et al. 1987; Groenewegen & Lamers 1989). Fluxes are normalized 
to the local continuum. The wind terminal velocity derived from the ob- 
served profiles is 3200 km s in both cases. The laboratory wavelengths of 
the doublet are indicated as vertical dotted lines. The figure indicates also the 
position of the violet shifted absorption components from the principal sys- 
tem at 1725 and 2850 km s~' in V 1974 Cygni ICassatella et al. 120041) . Such 
components cannot be clearly distinguished in GQ Mus due to the poor signal 
quality (the spectrum is the result of a 1 1-point running smoothing average). 
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FIG. 5. — Schematic representation of the energy spectrum of GQ Mus 
according to our model. The passbands of the photometric filters used in 
this work is indicated. The UV 1455A and supersoft X-ray (0.1—2.4 keV) 
fluxes are calculated from a blackbody spectrum, while the V, J, H, and K 
magnitudes are from a free-free emission spectrum as shown in the figure. 
Note that the flux scale of the dereddened blackbody and of the free-free 
emission corresponds to our 0.7 Mq white dwarf model on day 1 80 as shown 
in Fig|7] assuming a reddening of E(B — V) = 0.55. Horizontal short lines 
denote the energy flux in each filter on day 180; some of them are interpolated 
from other observational points. The soft X-ray flux is negligibly small on 
that day. The V and J magnitudes are larger than those of free-free because 
these bands are heavily contaminated by strong emission lines as explained 
in the Appendix. 

assumed that the spectral slope remained the same in the two 
spectra. The value so obtained, reported as an open circle in 
Figure [2] likely represents an upper limit to the true value be- 
cause, as shown in the same figure, the UV color was compar- 
atively cooler in the first observations. Figure [2] reports also, 
for comparison, the visual light curve obtained from the Fine 
Error Sensor (FES) counts on board IUE, once corrected for 
FES time dependent sensitivity degradation (see Cassatella et 
al. 2004 for details on the FES calibration). 

It appears from Figure [2] that the gradual fading of the 
UV and visual fluxes was interrupted, around JD 2,445,499 
(day 151), b y a se condary outburst, as also noticed by 
lHassall et al. I 0990). The secondary outburst had actually 
the appearance of a "UV flash" because of its especially large 
amplitude at short wavelengths. Indeed, compared with the 
IUE low resolution observations obtained just before and af- 
ter this event (days 108 and 202), the UV flux increased by a 
factor of 9 at 1455 A and by a factor 2.2 at 2885 A, while the 
visual flux increased only by a factor of 1.5. The consequent 
hardening of the UV continuum is also reflected by the up- 
turn of the C(1455-2885) color index, as shown in the same 
figure. 

3.3. Evolution of the UV emission lines 

Also the UV emission lines suffered from important 
changes in coincidence with the secondary outburst of day 
151, as shown in Figure [3] which reports, as a function of 
time, the flux in the most prominent permitted UV emission 
lines, as measured by us from the available IUE low resolu- 
tion spectra. The figure shows that, on day 151, the flux of 
the high ionization emission lines (He II 1640 A, N V 1240 A 
and C IV 1550 A) had decreased substantially while, on the 
contrary, it increased in the low ionization resonance lines of 
C II and O I, so implying a drastic change of the ionization 
structure of the emitting regions. 

The other important feature of the secondary outburst was 
the appearance of a strong P Cygni profile in the Si IV dou- 
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TABLE 1 
Chemical Abundance by Weight 



object H CNO Ne Na-Fe reference 

Solar composition 0.7068 0.0140 0.0018 0.0034 Grevesse & Anders (1989) 

GQ Mus 1983 0.37 0.24 0.0023 0.0039 Morisset & Pequignot (1996a 

GQMusl983 0.27 0.40 0.0034 0.023 Hassall et al. (19901 

GQMusl983 0.43 0.19 ■■■ ■■■ Andreae & Drechsel (1990) 



blet AA 1393.74, 1402.77 A, shown in Figure H The fig- 
ure reports also, for comparison, the Si IV profile of the 
Neon Nova VI 974 Cyg obtained 45 days after the outburst. 
A comparison of observed with theo retical P Cygni pro- 
files computed with the SEI method (lLamers et al. 1 119871: 
iGroenewegen & Lamers l!T989t) indicate that the terminal ve- 
locity of the wind was about the same for the two novae (« 
3200 km s -1 ). The special interest of the observations in Fig- 
ure [4] is that GQ Mus had entered the nebular phase well be- 
fore day 151, so that the presence of P Cygni profiles was 
not expected to be detectable at this stage. Indeed, P Cygni 
profiles were not detectable in an earlier IUE high resolution 
spectrum of day 49, despite the factor of two longer expo- 
sure time, nor they were at later stages. This strongly sug- 
gests that a short duration mass ejection episode took place 
around day 151. This episode caused efficient recombination 
to take place in the emitting region, observed as a substantial, 
although transitory, fainting of the high ionization emission 
lines (cf. Figj3). The fast recovery of the emission line spec- 
trum by day 202 suggests that the ejected mass must have 
been small. Also the presence of a flare at UV wavelengths 
rather than in the optical (cf. Fig. |2]i requires a small opacity 
in the UV and then a small ejected mass. 

4. MODELING OF NOVA OUTBURSTS 

We present a unified model for the IR, optical, UV, and su- 
persoft X-ray light curves of the 1983 outburst of GQ Mus. 
As in Paper I, our models are based on t he optically thick 
wind t heory of nova outbursts described in iKato & Hachisu I 
(fl994l) . 

4.1. Optically thick wind model 

After a thermonuclear runaway sets in on a mass-accreting 
WD, its envelope expands greatly to > 100 Rq, where 
R p h is the photospheric radius, and it then settles onto steady- 
state regime. The decay phase of the nova can be followed 
by a s equence of steady state solutions (e.g jKato & Ha chisu 
1994). Using the same m ethod and numerical techniques as in 
Kato& Hachisu (1994), we have followed the nova evolution 
by connecting steady state solutions along the envelope mass- 
decreasing sequence. 

The equations of motion, radiative diffusion, and conserva- 
tion of energy are solved from the bottom of the hydrogen-rich 
envelope through the photosphere under the condition that the 
solution goes through a critical point of a steady-state wind. 
The winds are accelerated deep inside the photosphere so that 
they are called " optically thick winds." W e have used updated 
OPAL opacities (Iglesias & Rogers Tl996l) . One of the bound- 
ary conditions for our numerical code consists in assuming 
that photons are emitted at the photosphere as a blackbody 
with a photospheric temperature T p h. X-ray and UV fluxes are 
estimated directly from the blackbody emission, but infrared 
and optical fluxes are calculated from free-free emission by 
using physical values of our wind solutions as mentioned be- 
low in 34.31 Physical properties of these wind solutions have 



already been published (e.g IHachisu & Kato 11200131 12004 ; 
I Hachisu et aU[l996l [l9993"bl 120001 l2003t IKato II 1 9831 [1997 . 
Il999l) . 

Optically thick winds stop after a large part of the envelope 
is blown in the winds. The envelope settles into hydrostatic 
equilibrium, while its mass decreases in time by nuclear burn- 
ing. Then we solve the equation of static balance instead of 
the equation of motion. When the nuclear burning decays, the 
WD enters a cooling phase, in which the luminosity is sup- 
plied with heat flow from the ash of hydrogen burning. 

4.2. Multiwavelength light curves 

In the optically thick wind model, a large part of the en- 
velo pe is ejected continuous ly during a relatively long period 
(e.g jKato & Hachisu 1119941) . After maximum expansion, the 
photospheric radius gradually decreases keeping the total lu- 
minosity (L p h) almost constant. The photospheric temperature 
(T p h) increases with time because of L pn = 4irRp h aTp h . The 
maximum emission shifts from the optical to supersoft X-ray 
through ultraviolet (UV) and extreme ultraviolet (EUV). This 
causes the luminosity to decrease in the optical and to increase 
in the UV, until it reaches a maximum. The following decay in 
the UV is accompanied by an increase of the supersoft X-ray 
range. These timescales depend crucially on the WD parame- 
ters such as th e WD mass a nd the chemical composition of the 
envelope (e.g. lKatolll997l Paper I). Thus, we can follow the 
development of optical, UV, and supersoft X-ray light curves 
by a single modeled sequence of steady wind solutions. 

4.3. Free-free light curves 

Spectra of novae show blackbody features at very early 
stages, but free-free emission fr om optically thin plasma 
will eventually dom inate (e.g., Gallagh er & Nev I 1 19761 : 
iRrautter et al. Ifl9 84l). During the optically thick wind phase 
(see Paper I for more details), the extended regions outside the 
photosphere are optically thin. The optical and IR free-free 
emission fluxes arising from these regions can be estimated 
from 

F v oc / N e NidV oc / -^^r 2 dr oc (2) 

J JRph V willed V ph P h 

where F v is the flux at the frequency v, N e and A/, are the 
number densities of electrons and ions, respectively, V is 
the volume of the optically thin region, /? pn is the photo- 
spheric radius, M w j n d is the wind mass loss rate, and v pn is 
the photospheric velocity. In equation © we assume that 
the electron temperature and the degree of ionization are con- 
stant in the free-free emitting region. We also assume that 
N e oc p„i„d, Ni oc p w i n d, and use the continuity equation, i.e., 
Pwind =Mwind/47rr 2 v wind , where p wind and v wind are the density 
and velocity of the wind, respectively. Finally, we assume that 
Vwind = const. = v pn in the optically thin region. 

After the optically thick wind stops, the total mass of the 
ejecta remains constant in time. The flux from such homolo- 
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log[days after outburst] 

FIG. 6. — Calculated light curves of free-free V magnitude, UV 1455 A, and supersoft X-ray (0.1—2.4 keV) fluxes are compared with the observations. Open 
diamonds: ROSAT X-ray (0.1-2.4 keV) observations from iShanley et al. 1 11993) and Orio et al. (2001). Open circles: WE UV 1455 A ob servations. Open 
triang les: visual and V magnitudes are from IAU Circulars 3764, 3766, 3771, 3782, and 3853. Open squares: V magnitude observations from Whitelock et al. 
1 1984) and IUE V*pes data. Small circles: visual magnitude observations from RASNZ and AAVSO. Our best-fit model consists of the white dwarf with mass 
0.75 Mq for the envelope chemical composition of X = 0.55, Xcno = 0.20, and Z = 0.02. An arrow indicates the epoch when the optically thick wind stops. The 
upper thin solid line is the same as the thick solid line, but lifted up by 4.5 mag to represent contribution of strong emission lines to the V bandpass (Papers I and 
II). Model predictions for white dwarf masses Mwd = 0.7 Mq (thin solid line) and Mwd =0.8 Mq (dash-dotted line) are also added for comparison. 
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FIG . 7 . — Same as in Fig. [6] but for our best-fit model of Mwd =0.7 Mq for the envelope chemical composition of X = 0.45, Xcno = 0.35, and Z = 
other cases of the white dwarf mass, Mwd =0.65 Mq (thin solid line) and Mwd =0.75 Mq (dash-dotted line) are also added for comparison. 
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gously expanding ejecta is approximately given by 



N e NidV cx p 2 ej V ej cx ( ^ 

Vej 



" 3 ocr 3 , 



(3) 

where p e j, V e j, and M e j(= const.) are the density, volume, and 
total mass, respectively, of the ejecta. We assume that the 
ejecta are expanding at a constant velocity, v. So we have 
the radius of the ejecta of R = vt, where t is the time after the 
outburst. The proportionality constants in equations (O and 



(01 cannot be determined a priori because radiative transfer is 
not calculated outside the photosphere: these were determined 
using the procedure described below in §|5] 

4.4. System parameters of optically thick wind model 

The light curves of our optically thick wind model are pa- 
rameterized by the WD mass (M WD ), the chemical composi- 
tion of the envelope (X, Xcno, X Ne , and Z), and the envelope 
mass (AM enV; o) at the time of the outburst (JD 2,445,348.0). 
For the metal abundance we adopt Z = 0.02, which also in- 
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FIG. 8. — Same as Fig. [6] but for our best-fit model of the white dwarf mass of 0.65 Mq for the envelope chemical composition of X = 0.35, Xcno = 0.30, and 
Z = 0.02. Two other cases of the white dwarf mass, Mwd =0.6 Mq (thin solid line) and Mwd =0.7 Mq (dash-dotted line) are also added for comparison. 
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FIG. 9. — Near-infrared multiwavelength light curves (/, J, H, and K) for free-free emission, based on equation (2)> and our model light curve corresponding 
to Mwd = 0.70 Mq, X = 0.45, and Xcno = 0.35, together with the observations. Thick solid line: fit of the early phase. Thin solid line: fit of the late phase 
after the transition phase (see, e.g., Papers I and II). Open circles: end of the wind phase, (a) / magnitudes from free-free emission: These th ick and thin lines 
are overlapped, (b) J magn itudes of free-free emission. The He I A10830 contributes to the J band from the very early phase of the outburst i Whitel ock et al. I 
H984t lKrautter et al. 1984), so one more line (dash-dotted) is added to rep resent the very early pha se, ( c) H magnitudes of fre e-free emission, (d) K magnitudes 
of free-free emission. Observational data (/, J, H, and K) are taken from Whitelock et al. 1 1984) and Krautter et al. 1 1984). Essentially the same figures are 
obtained for the other two cases of M W d = 0.75 M , X = 0.55, X CN0 = 0.20 (see Fig. |6j, and M WD = 0.65 Mq, X = 0.35, and Z CNO = 0.30 (see Fig. [8). 
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eludes carbon, nitrogen, oxygen, and neon with solar compo- 
sition ratios. 

Three different sets of abundance determinations for the 
ejecta of GQ Mus are available from the literature, but their 
values are rather scattered, as it appears from TableQ] We here 
consider three cases, in order of decreasing hydrogen content 
in the envelope, not directly corresponding to the three sets of 
abundance determination in TableQ] a) X = 0.55, Xqno = 0.20, 
X Ne = 0.0, and Z = 0.02; b) X = 0.45, X CN0 = 0.35, X Ne = 0.0, 
and Z = 0.02; c) X = 0.35, X CN o = 0.30, X Ne = 0.0, and 
Z = 0.02. These composition sets correspond to assuming 
25%, 55%, and 100% mixing of C+O WD matter with a 
hydrogen-rich envelope with solar composition. If we change 
^cno while keeping the hydrogen content X constant, the light 
curves hardly change provided Xcno ^ 0.2. 

We have searched for the best fit model by changing the 
WD mass in steps of 0.05 Mq for the above three sets of 
chemical compositions. 

5. LIGHT CURVE FITTING 

We apply to GQ Muscae 1983 the model light curves of 
classical novae described in the previous section , and evaluate 
its fundamental parameters. The optical V magnitude and the 
near infrared /, J, H, and K magnitudes are calculated from 
free-free emission, while the UV 1455 A and supersoft X-ray 
fluxes are obtained from blackbody emission, as illustrated in 
Figure [5] 

Models and observations are compared in Figures |6] — [8] 
which refer to, respectively, the three following chemical 
compositions sets: (X,Y,X C no,Z) = (0.55,0.23,0.20,0.02), 
(0.45, 0.18, 0.35, 0.02),and (0.35,0.33,0.30,0.02). The white 
dwarf mass range explored is Mwd = 0.60 to 0.80 Mq, in 
0.05 Mq steps. In each figure, the thick solid line represents 
the white dwarf mass value that leads to the best representa- 
tion of the data for the above chemical composition sets. 

5.1. Supersoft X-ray fluxes 

The decay time of the supersoft X-ray flux is a good in- 
dicator of the WD mass (Paper I). Through a careful vi- 
sual inspection to the data, we have s elected the models that 
best fitted the ROSAT observations (IGgelman et al. 1 119931: 
Shanlev et al~l 11995b lOrio et al. I [20011: Barman & Krautterl 
20011) . The model and the observed supersoft X-ray light 
curves are reported in Figures [6]- [8] 

5.2. Continuum UV 1455 A fluxes 

The light curves in the UV 1455 A continuum are a good 
indicator of photospheric temperature during the early decay 
phase of no vae. Such curves are in general quite smooth, 
as shown by Cassat elTa et al. I (120021) for several objects ob- 
served with IUE. A comparison of models with observed UV 
light curves (see Figs. |6]-[8]l shows that this is also the case 
of GQ Mus, except for the two earliest observations on day 
37 and 49, and that of day 151. Leaving these data apart, we 
find that the model that best fits the data is that of a white 
dwarf with mass M WD = 0.75, 0.7, and 0.65 Mq, for the three 
above chemical composition sets, respectively. As for the ob- 
servation of day 151, the large deviation from the model's 
smooth trend is clearly due to the occurrence of the UV flash 
discussed in §3.2. The otherwise good agreement between 
models and observations shown in the above figures may sug- 
gest that also the observations of day 37 and 49 are due to a 
UV flash (see discussion in j 
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FIG. 10. — Evolution of wind mass loss rate (M w i n( j: dashed line), pho- 
tospheric temperature (7^: upper thin solid line), photospheric radius (Rph : 
middle thin solid line), photospheric luminosity (L p i,: lower thin solid line), 
and X-ray flux (thick solid line) of our Mwd = 0.7 Mq model for the enve- 
lope chemical composition of X = 0.45, ^cno = 0-3 5, and Z = 0.02. Open dia- 
monds: Observational X-ray count rates taken from Sha nlev et al. 1 (1995 ) and 
IQrio et al. I f200ll) . Open squares and open triangles: photosp heric temper- 
atures and photospheric radii, respectively, taken from Balman & Krautter 
12001). 



5.3. Optical and infrared fluxes 

The visual light curves computed from free-free emission 
are compared with the observed ones in Figures [6]- [8] for the 
above chemical composition sets. As discussed in §1 and in 
Papers I and II, visual magnitudes are contaminated by strong 
emission lines which will eventually dominate over the con- 
tinuum, causing an increasing deviation from our free-free 
models. In GQ Mus, the forbidden [O III] AA 4959, 5007 
emission fines alread y appeared 39 days after the outburst 
dKrautter et al. | [l984). At about this date the observed visual 
light curve did actually start to show an increasing deviation 
from the free-free emission model until, by about day 500, the 
contribution from the emission lines stabilized so that the ob- 
served light curve recovered the the shape of the model curve. 
This effect is approximatively taken into account by lifting 
up the template light-curve of GQ Mus by about 4.5 mag, as 
shown in Figures|6]-|8] 

On the other hand, the near infrared IJHK bands are not so 
heavily contaminated by emission lines, as shown in Figure 
|9l whe re the IJHK photometric data from Whi telock et al. I 
(1984) are compared with our theoretical light curve. The 
figure shows that the model fits the observations reasonably 
well until about day 100 (thick solid line in the figure), being 
anyhow the deviations after that date sensibly smaller than in 
the visual. The contribution from emission lines is particu- 
larly small in the /, H, and K bands, whereas the emission 
line of He I A10830 c ontributes somewhat to t he J band from 
the v ery early phase (Whitelock et al. 19844 iKrautter et al. I 
1984). To mimic the different contribution by the emission 
lines at different phases of the nova development, it is suffi- 
cient to lift up the template free-free line two times (starting 
from the dash-dotted line representing the very early phases, 
to the thick solid line representing the intermediate phases, 
and then finally to the thin solid line, representing the late 
phases). 

To summarize the results in this section, the white dwarf 
mass that best reproduces the X-ray, UV, optical, and near IR 
is Mwd = 0.7 ± 0.05 Mq, where the error bar reflects mainly 
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FIG. 11. — Photospheric temperature (T p t,: upper thin solid line), pho- 
tospheric luminosity (L p i,: middle thin solid line), photospheric radius (R p i,: 
thin solid line from left-upper to right-lower), and UV 1455 A flux (thick 
solid line) of our Mwd = 0-7 Mq model for the envelope chemical composi- 
tion of X = 0.45 and X CN o = 0.35. Open circles: UV 1455 A fluxes. Visual 
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FIG. 12. — Visual and V magnitudes during 300 days after the optical max- 
imum. The same symbols as in Figs. |6]-[8]except the WE V FES magnitude 
(filled squares). Two free-free light curves are plotted both for the IUE Vfes 
magnitudes and for Whitelock et al.'s (1984) V magnitudes. The Vfes data 
are 0.5-0.8 mag brighter than those of Whitelock et al.'s. The ?3 time (left ar- 
row) is estimated to be 50 days from Whitel ock et al.'s V m agnitudes together 
with the peak brightness of 7.2 mag i Liller & Overbeek 1983). Note that the 
decay time (right arrow: ?3 = 122 days) corresponding to our universal de- 
cline law (lower solid free-free line) is sensibly larger than that previously 
quoted. See text for more details. 

the uncertainty on chemical composition. 

6. DISCUSSION 

6.1. Overall Development of the Nova Outburst 

In this subsection, we analyze the overall development of 
the nova outburst based on our model of Mwd = 0.7 Mq , X = 
0.45, Xcno = 0.35, and Z = 0.02 (Fig. Figures M - E] 
report some relevant model predictions and observations as a 
function of time. First we point out two important epochs of 
the nova development: one is the end of optically thick wind 
phase on day 1000, and the other is the end of hydrogen shell 
burning on day 3300 as summarized in Table [3 Because of 
the rapid shrinking of the the photospheric radius (K p h), also 
the photospheric temperature (T p h) rapidly increases until the 
wind stops, as indicated in Figure[l0] Also the wind mass loss 
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FIG. 13. — Color-temperature relation for nova wind solutions. The UV 
color of C(1455-2885) is defined by -2.51og(F(1455)/F(2885)). Filled 
circles connected with a solid line: reddening-corrected C(1455 — 2885) 
v.s. the photospheric temperature T c ff corresponding to our best fit model. 
Open triangles: color- temperature relati ons estimated from Fig. 10 of 
IHauschildt etaTI (199% and Fig.8 of IShort et al~l (199% . See text for more 
details. 
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FIG. 14. — Distance-reddening relations derived from the UV 1455 A 
fitting (solid line labeled "UV 1455 A ") of our model M W d = 0.7 Mq, 
X = 0.45, and Xcno = 0.35 and the maximum magnitude vs. rate of decline 
(labeled "MMRD1" and "MMRD2"). The dashed line labeled "MMRD1" is 
a MMRD relation calculated from Schmidt-Kaler's law, i.e., eq. ^5} together 
with ?3 = 50 days. The solid line labeled "MMRD2" is the same MMRD rela- 
tion calculated from /3 = 122 days estimated from our universal decline law. 
An arrow indicates our value of E(B-V) = 0.55, which is an error-weighted 
mean of E(B-V ) = 0.51 ±0.06 and 0.58 ±0.04 in §3.1. See text for more 
details. 

rate (M w ; nc j) decreases rapidly from ~2x 10~ 4 M© yr -1 to 1 x 
10~ 7 M Q yr" 1 during the first 1000 days. Thus the total mass 
ejected by the winds amounts to ~ 2 x 10~ 5 M Q . After the 
wind stops, the photospheric radius becomes smaller than ~ 
0. 1 Rq, and the photospheric temperature eventually becomes 
larger than 20 eV, so that the nova enters a supersoft X-ray 
phase. After the hydrogen shell burning ends, the nova cools 
down rapidly, and a fast decay phase of the supersoft X-ray 
flux follows, as already shown in the previous section. In what 
follows, we summarize some relevant observational aspects 
and make a comparison with our model predictions. 
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From the coronal lines in the optical spectrum Diaz et al. 
(1 19921) deduced a luminosity of 10 37 - 10 38 ergs s" 1 and a tem- 
perature of (2-3) x 10 5 K. These values are consistent with 
our model of a 0.7 M© WD, which provides r p h = 4 x 10 5 K 
and Lph = 6 x 10 37 ergs s" 1 about 2700 days after the outburst, 
a s it may be apprecia ted from Figure [TOl 

IDiaz et aL~l d 19951) estimated a central source temperature 
of 164,000 K, a luminosity of 350 Lq, and an ejecta mass of 
5 x lO^M©, about 4000 days after the outburst. For our best- 
fit model of a 0.7 M© white dwarf with an initial envelope 
mass of AMo = 2.7 x 10~ 5 M©, the total mass lost by the op- 
tically thick wind is AM wind = 1.9 x 10~ 5 M©. We also obtain 
7p h = 1.5 x 10 5 K and L ph = 4 x 10 35 ergs s~ Y ks 100 L©, 4000 
days after the outburst, as shown in Figure [10] Our values 
are the n very consistent with those obtained by IDiaz et aT71 
(1 19951) . except for their larger ejected mass (see below for an- 
other estimate). 

The envelope mass itself depends mainly on the white 
dwarf mass and slightly on the chemical composition. For 
a 0.65 M© white dwarf, X = 0.35, and X C no = 0.30, we ob- 
tain a slightly larger envelope mass, AM,, = 3.6 x 10 _5 M Q , 
and a larger wind mass loss, AM w i nc j = 2.3 x 10~ 5 M©. On 
the other hand, for a 0.75 M© white dwarf with X = 0.55, 
and Xcno = 0.20, we obtain a slightly smaller envelope 
mass, AMo = 2.6 x 10~ 5 M©, and a smaller wind mass loss, 
AMwind =1.7 x 10- 5 M©. 

Usin g their photoionization model, Morisset & Pequignot 
(1996b) estimated the temperatures during the static phase 
and at the hydrogen burning turnoff to be 2.5 x 10 5 K and 
4.1 x 10 5 K, respectively, and an ejected shell mass of 8 x 
10~ 5 M©. These values are consistent with our model, i.e., 
> 2.4 x 10 5 K, > 4.4 x 10 5 K, and 2 x 1O _5 M . The above 
authors determined the duration of the wind phase and of the 
static hydrogen burning phase to be < 0.52 yr and 8.8 yr, to 
be compared with our values of 2.7 yr and 6.3 yr, respectively. 
The much longer duration of the wind phase obtained by us 
is likely due to the rapid decrease of the wind mass loss rate 
with an e-folding time of 0.24 yr, corresponding to a decrease 
by a factor of 10 in 0.7 yr. 

6.2. Photospheric Temperature Development 

Figure [13] shows the reddening corrected UV color index 
C(1455-2885) as a function of the photospheric temperature 
obtained from our best-fit model at the time of the individual 
observations (filled circles). The figure indicates that the color 
index decreases smoothly with increasing temperature except 
for its abrupt decrease at the time of the UV flash on day 151 
(note that the temperature increases with time in these early 
phases). 

It is interesting to compare the co lor- 7ph values reported in 
Fig ure [131 with t hose derived from Hausch ifdt et al. I (1 19971) 
and S hort et al. I (1 19991) . who calculated early-time nova spec- 
tra using non-LTE atmosphere codes with winds having den- 
sity and velocity laws within the nova envelope. For this com- 
parison, we have taken the values o f C(1455-2885) from Fig- 
ure 10 of lHauschildt et al. I (1 19971) for T eS = 15, 000, 20,000 
25,000 , and 30,000 K, and from Figure 8 of IShort et al. I 
(1999) for T eff = 35,000 K. The model values of C(1455- 
2885) so obtained are plotted as a function of the correspond- 
ing effective temperature in Figure Qj] (open triangles). It 
clearly appears from Figure Q~3] that the two sets of values 
are very consistent with each other, except for the UV flash 
on day 151. This result supports our assumption that the UV 



1455 A flux is reasonably well accounted for by blackbody 
emission. Indeed, line blanketing is rather small at this wave- 
length, so that the blackbody model does not deviate much 
from Hauschildt et al.'s non-LTE models, as already discussed 
in Paper I. 

6.3. Emergence of the secondary component 

The mass of the donor star (the sec ondary component) can 
be es timated from the orbital period. Diaz & Steinerl (1 19891 
119941) obtained P orb = 0.05936 days (1 .425 hr) from the orbital 
modulations with an amplitude of 0.6 mag. Using Warner's 
(1995) empirical formula 

^ « 0.065 f-^V /4 , forl.3<-^<9 (4) 
M© \ hours / hours 

we get M2 = 0.10 M©. The orbital separation is then a = 
0.59 Rq for M W d = 0.7 M©, the effective radius of the 
Roche lobe for the primary component (white dwarf) is 
R\ = 0.33 Rq, and the effective radius of the secondary is 
R2 = 0.14 Rq. In our model, the companion emerges from 
the white dwarf envelope when the photospheric radius of 
the white dwarf shrinks to /? p h ~ 0.6 Rq (the separation) or 
0.3 Rq (the Roche lobe). This happens on about day 330 or 
500, respectively, in our best-fit model with Mwd = 0.7 M©, 
as shown in Fi gures [TOl and ITTI 

According to White lock et al. I (119841) . strong infrared coro- 
nal lines appeared 57 days after optical maximum in VI 500 
Cyg, but were not present in GQ Mus as late as 97 days af- 
ter maximum. Their appearance roughly coincides with the 
emergence of the companion from the white dwarf photo- 
sphere: this suggests that these strong coronal lines arise from 
the shock between the white dwarf wind and companion star. 
Indeed, the emergence of the companion in VI 500 Cyg took 
place about 50 days after optical maximum (Paper I), which is 
consistent with the appearance of the line 57 days after max- 
imum. In the case of GQ Mus, the absence of strong coronal 
lines on day 97 days is consistent with our model in which the 
companion emerged about 33 or 500 days after maximu m. 

In their study of GQ Mus. iKrautter & Williamsl (1 19891) re- 
ported that the [Fe X] A6374 coronal line first appeared 2.2 
to 3.4 yrs after the outburst, and became the strongest about 
4 years after it. They argued that the coronal line was due 
to photoionization from a hot radiation source rather than to 
collisional excitation. In our 0.7 M© white dwarf model, the 
optically thick wind stopped about 960 days (2.7 yrs) after the 
outburst and then the photospheric temperature gradually in- 
creased to above 20 - 25 e V 4 yrs after the outburst, as shown 
in Figures[lO]and[TT] These authors suggested that ~3 x 10 5 
K in 1984, which is consistent with our photospheric temper- 
ature of 7p h = 2.2 x 10 5 K 2.2 yr after the outburst (Figs. [10] 
andQj}- 

6.4. Hard X-ray component 

Ogelman et al. I (119871) reported that the count rate of the 
EXOSAT low-energy telescope was about 3 x 10" 3 c s" 1 400 
days to 700 days after optical maximum, a value that gradu- 
ally decreased to 1 x 10~ 3 c s" 1 about 900 days after it. 

If the X-rays originates from the shock between the wind 
and the companion star, the soft X-ray behavior is consistent 
with our 0.7 M© white dwarf model (as well as 0.75 M© and 
0.65 Mq white dwarf models), because the emergence of the 
companion star is predicted to happen on day 330, and the 
wind stopped at 960 days after optical maximum. 
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TABLE 2 
Properties of GQ Muscae 



subject 


data 


units 


reference 


discovery 


. 2,445,352.6 


JD 


Liller & Overbeek (1983) 


nova speed class 


moderately fast 




Whitelock et al. (,1984) 


h 


. 18 


days 


Whitelock et al. (1984) 


h 


. 48 


days 


Whitelock et al. (1984) 


Mv.peak from r 3 


. -7.55 


mag 


Whitelock et al. (1984) 


distance from 


. 5 


kpc 


Whitelock et al. (1984) 


dust 


no 




Whitelock et al. (1984) 


orbital period 


. 1.425 


hr 


Diazetal. (1995) 


E(B-V) 


. 4.5 




Krautter et al. (1984) 


wind phase 


. < 0.5 


yr 


Morisset & Perjuienot (1996b) 


H-burning phase 


. 3300 


days 


Shanleyetal. (1995) 



TABLE 3 
Summary of our model 



subject 


data 


units 


outburst day 


. 2,445,348.0 


JD 


opt. maximum 


. 2,445,352.6 


JD 


h 


. 67 a 


days 


h 


. 122 a 


days 


M v . pcilk from f 3 


. -6.53 b 


mag 


distance from % 


. 5.0±0.3 


kpc 


secondary mass 


. 0.1° 


M 


E(B-V) 


. 0.55 ±0.05 




distance by UV fit 


. 4.9 ±0.9 


kpc 


^break 


. 275 


day 


WD mass 


. 0.7 ±0.05 


Mq 


WD envelope mass 


. (2.6-3.6) x 10- 5 


Mq 


mass lost by wind 


. (1.7-2.3) x 10- 5 


Mq 


wind phase 


. 1000 


days 


H-burning phase 


. 3300 


days 


separation 


. 0.6 


Re 


companion's emergence 


. 330 ±30 


days 


UV flash luminosity 


. >45,000 e 





a 12 and ?3 are calculated from our fitted universal decline law. cal- 
culated from equation (5) together with ti above. c estimated from 
equation see Paper I. e estimated from a blackbody of T ~ 

10 5 K. 

6.5. Distance 

IWhitelock et al. I (11984 and lKrautter et al. I (11984 have es- 
timated the distance to GQ Mus from the apparent and the ab- 
solute V magnitudes at maximum, my and M v . This latter was 
computed from the tj, time through the Schmidt-Ka ler Maxi- 
mum Magnitude Rate of Decline (MMRD) relation (ISchmidt I 
fl957h : 

M v =-11.75 + 2.51ogf 3 . (5) 

Despite the sensibly different input values adopted of 48 and 
40 days for f 3 , 0.4 and 0.45 for E(B-V), and 7.2 and 6 for m v , 
these authors obtained quite similar values for the distance, 5 
kpc and 4.8 kpc, respectively. 

A re-analysis of the merged visual magnitude data present 
in the literature (see p lot in Fig. fT2l confirms ft ~ 50 days, in 
close agreement with Whitel ock et al. I (11984. so that My 



7.5 [ cf. equation 10]. If, together with IWhitelock et al. I 
( 1984), we take mymsa. = 7-2 as the best estimate of the ap- 
paren t magnitude at optical maximum dLiller & Overbeek] 
1983), we finally obtain the following relation between dis- 
tance and E(B-V): 

m v -M v = -5 + 51og d + 3AE(B-V)= 14.7, (6) 

which is labeled "MMRD1" in Figure [14] In particular, the 
value E(B-V) = 0.55 ±0.05 obtained in §3.1, corresponds to 
a distance of d = 4.0 ± 0.3 kpc. 
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FIG. 15.— Energy distribution for three bands of UV 1455 A, 2885 A, and 
IUEFES visual (open triangles) on day 151. Blackbody energy distributions 
are also plotted for three temperatures of 50,000, 100,000, and 200,000 K. 

Following the same proced ure outlined in Paper I and 
iKato & Hachisu I (l2005Ll2~007l) . an independent estimate of 
the distance to GQ Mus can be obtained by comparing the 
observed light curve in the 1455 A continuum with the cor- 
responding model fluxes in Figures [6] - [8] HH and [T2] The 
calculated flux at A = 1455 A at a distance of 10 kpc for our 
adopted model (0.7 M Q ,X = 0.45, Xqno = 0.35) is F™ od = 2.85 
xlO" 12 ergs cm" 2 s" 1 A" 1 on JD 2,445,455.6. The observed 



flux at the same date is F£ hs = 1.78 x 10 13 ergs cm 2 s 



A" 



From these values one obtains the following relation between 
distance and reddening: 

d 



^obs 



Mr d = 51og(- 



-)+A x E(B-V) = 3.01, (7) 



10 kpc 

where m A = -2.51og(f\), and A\ = 8.3 dSeaton 1119791) . This 
curve, labeled "UV 1455 A," is reported in Figure[l4] In par- 
ticular, equation (|7]i provides a distance of 4.9 ± 0.9 kpc for 
E(B-V) = 0.55 ± 0.05, which is sensibly larger than that ob- 
tained from MMRD 1 . 

An alternative way to determine the distance of GQ Mus 
is to compare the shape of the theoretical visual light curve 
(based on the free-free model) with the observed visual 
fluxes. Such a comparison, done in Figure [7] is very instruc- 
tive because it shows that the object was up to 1.5 magnitudes 
brighter than predicted by our universal light curve model dur- 
ing the earliest phases (until day 10). It is then legitimate to 
consider GQ Mus as a super-bright nova. An extreme case of 
this type was VI 500 Cyg, which, near peak luminosity was 
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about 4 magnitudes brighter than the E ddington limit for a 
1.0 M Q white dwarf (iFerland et al. 111986b . It is interesting to 
note that if such a super-bright phase is ignored in GQ Mus, 
the ?3 and ?2 times of the theoretical light curve that best fits 
the observations are considerably longer than the observation 
indicates at face value: we in fact estimate £3 = 122 days and 
?2= 67 days. From equation (0 we then obtain Mff max = —6. 53 
at the time of maximum light (t = 8 days). Since the apparent 
magnitude at visual maximum of our theoretical light curve is 
Mffmax = 8.66 mag (see Fig. |7), the distance to GQ Mus can 
be estimated to be 5.0 ± 0.3 kpc for E(B-V) = 0.55 ±0.05. If 
the reddening is not fixed, the following relation applies: 

m ff ,max-M ff , max =-5 + 5\ogd + 3AE(B-V) = 15.19, (8) 

which is plotted in Figure [14] (labeled "MMRD2"). The two 
curves UV 1455 A and MMRD1 cross at d = 3.5 kpc and 
E(B-V) = 0.64, whereas the curves UV 1455 A and MMRD2 
cross at d = 5.1 kpc and E(B-V) = 0.54, as it appears from 
Figure Since E(B-V) = 0.55 is the weighted mean of 
the reddening determinations done in §3. 1, we take this value 
together with d= 5 kpc as the best estimates compatible with 
the observations. The errors on the reddening and the distance 
are ~ 0.05 dex and ss 0.5 kpc, respectively. 

6.6. The UV flashes 

Here we estimate the total luminosity of the UV flash on 
day 151 and discuss its outburst nature. From Figure 12 we 
have obtained the excess of energy above the smooth de- 
cline for each of the three bands, i.e., 1455 A, 2885 A, and 
Vfes, which are plotted in Figure [15] These fluxes are con- 
sistent with the energy distrib ution from a blackbod y with a 
temperature of > 100,000 K. lHassall et aTl (1 19901) also es- 
timated the temperature from a Zanstra-like method, based 
on the He Il/H/3 ratio, to be 85,000 to 100,000 K. The pho- 
tospheric temperature of our 0.7 M Q white dwarf model are 
about 30,000 K (on day 108) before the UV flash and about 
60,000 K (on day 202) after the UV flash. 

Assuming blackbody emission, we can estimate the size 
of the emitting region from (R/df w 10" 21 for T = 100,000 
K. Taking d = 5 kpc, the radius of the emission region is 
R ~ 7 Rq. The total flux is estimated to be 1.7 x 10 38 erg s" 1 
(« 45,000 L Q ) from L = 4nR 2 aT 4 . On day 151, the photo- 
spheric radius of our 0.7 M Q white dwarf model is as small as 
1.3 Rq and the photospheric temperature is as low as 42,000 
K. Therefore, the estimated radius and blackbody tempera- 
ture, R ~ 7 Rq and T > 100,000 K, suggest an episodic ex- 
pansion and strong heating of the photosphere. This was prob- 
ably due to a dynamical mass ejection episode and to strong 



shock heating. The UV flash is anyhow a short-lived event 
that hardly contributes to the total ejecta mass, as suggested 
by the fast recovery of the emission line spectrum by day 202 
and by the negligible effect on the optical light curves. 

The first observation of IUE on day 37 also shows a high 
UV flux compared with our white dwarf model. A close 
look at the optical brightness in Figure Q~2] indicates a small 
peak around day 37. This event may also be due to another 
mass ejection episode. Note that these mass ejection episodes 
add to the major underlying continuous mass ejection pro- 
cess from the optically thick wind, which endures until day 
~ 1000, as shown in Figure [10] 

7. CONCLUSIONS 

We have applied the "universal decline law" of classical no- 
vae described in $4]to GQ Muscae 1983 and derived various 
parameters of the nova. Our main results may be summarized 
as follows: 

1 . We show that the "universal decline law" reproduces well 
the observed light curves of GQ Mus in the optical and in the 
near infrared /, J, H, and K bands. 

2. Our blackbody light curve model for the UV 1455 A 
band can reproduce the observed UV 1455 A fluxes except 
two UV flashes on day 37 and 151. 

3. The UV flash on day 151 described in jj3.2l was accom- 
panied by a mass ejection episode seen as a fast wind with a 
terminal velocity of about 3200 km s -1 . 

4. An analysis of the IUE reprocessed data of GQ Mus 
indicates E(B-V) = 0.55 ±0.05, a value that is larger than 
previously reported. 

5. We find that the mass of the WD component of GQ Mus 
is 0.7 ±0.05 Mq for the adopted envelope chemical composi- 
tion of X = 0.35 -0.55, Xcno = 0.2-0.35. This value has been 
derived by comparing predicted fluxes from our models with 
the observations in the supersoft X-ray, in the UV, and in the 
optical and near infrared, obtained at different times. 

6. We have estimated a mass of AM w j n d ~2x 10~ 5 M Q lost 
by the wind. 

7. We have estimated a distance of d ~ 5 kpc. 
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helped to improve the paper. This research has been supported 
in part by the Grant-in- Aid for Scientific Research ( 1 65402 1 1 , 
16540219, 20540227) of the Japan Society for the Promotion 
of Science and by INAF PRIN 2007. 



APPENDIX 

COLOR OF FREE-FREE LIGHT CURVES OF NOVAE 

In this and previous papers we have modeled the time evolution of novae in the continuum and compared the results with 
observed light curves. Such a study is hampered by the presence of emission lines, which give rise to sensible distortions in most 
photometric bands (see §1), especially at some nova phases. To evaluate the contamination by emission lines, we take advantage 
of the fact that the intrinsic colors of free-free emission are constant with time because of XF \ oc A -1 for free-free emission. Here, 
we consider four classical novae that are well observed in various photometric bands: V1500 Cyg, V1668 Cyg, V1974 Cyg, and 
GQ Mus. 

Intrinsic V-R, V-I, V-J, V-H, and V -K colors from free-free emission, as calculated from \F\ oc A -1 , are reported in the 
second column of Table[4] If reddening is known, observed colors are obtained from 

m v -m x = (M v -M x ) Q + c x E{B-V), (Al) 

where (M v -M\)q is the intrinsic color and c\ is the reddening coefficient listed in the second and third columns of Table [4] 
respectively. 




logfdays after outburst] 

FIG. 16. — Observed y, V, I, J, H, and K light curves for the classical nova VI 500 Cyg. Each light curve is shifted down to the y light curve by Am\ — AV = 
y — V = 0.3, AI = y— I = 2.2, AJ = y—J = 2.75, AH = y-H = 2.65, and AK = y—K = 3.1. These light curves are well overlapped between day 6 and 30 (days 
after outburst). Filled circles: y magnitudes. Asterisks: V magnitudes. Plusing open circles: / magnitudes. Filled triangles: J magnitudes. Open squares: H 
magnitudes. Open circles: K magnitudes. Solid line with a large open circle: our model light curve for the 1.15 Mq WD with X = 0.55, Xcno = 0. 10, Xn c = 0.03, 
and Z = 0.02 (see Paper I). The large open circle at the right-lower edge of the line indicates the end of optically thick wind phase. 
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In the following subsection, we will compare the observed colors with those calculated from equation (lAlb for classical 
novae V1500 Cyg, V1668 Cyg, V1974 Cyg, and GQ Mus. Figures [1611191 demonstrate that the shape of light curves is almost 
independent of the wavelength, which is a characteristic feature of free-free emission (see equation J2J). In these figures, we 
shift each observed light curve down by Am\ and overlap it on the observed y (or V) light curve. Then, we obtain the color of 
Am\ = m y - m\ (or A.m\ = mv-m\) from the definition of m\ + Am\ = m y (or m\ + Am\ = my). The colors thus obtained are 
listed in each figure and plotted in Figurel20l 
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V1974 Cyg (Nova Cyg 1992) 
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FIG. 18. — Same as Fig |16l but for the classical nova V1974 Cyg. No y and / magnitudes are available for this nova, so that we show only three colors of 
V —J, V —H, and V —K. Each light curve is overlapped with the V light curve during day 10-30. Solid line with a large open circle: our model light curve for the 
1.05 M WD with X = 0.55, X C no= 0.10, Z Nc = 0.03, and Z = 0.02 (see Paper I). Solid line withr 3 : a decline law of F x oc r 3 (seeeq.(3)). 
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FIG. 19. — Same as Fig |16l but for the classical nova GQ Mus. No y magnitudes are available for this nova, so that we show four colors of V -I, V—J, V -H, 
and V — K. Each light cur ve is overlapped with the V light curve during day 8-50 except for the J magnitudes. The / magnitudes are fitted only during day 
20-100. As mentioned in j|5.3l the J magnitude brightened up by 0.6 mag after day 15. If we overlap the J light curve between day 8 and 15, the color decreases 
from V— J = 3.1 to V — J = 2.5 as shown in the figure. Solid line with a large open circle: our model light curv for the 0.7 Mq WD with X = 0.45, Xcno = 0.35, 
and Z = 0.02 (see i|5}. Dashed line with a large open circle: same as that for the solid line but vertically shift down by 0.6 mag to match the J magnitudes during 
day 8-15. 



V1500 Cyg 

For VI 500 Cyg, we us e the y photometry in lLockwood & Millis I (1 19761) and the V data in [ empest 1 (fl979h. / d ata in 
The & van der Klisl (1 19761) . The J, H, and K observations were taken from lEnnis et al71 (1 19771) . Kawara et al. I (1 19761) . and 
Galla gher & Ney I (119761) . These light curves are plotted in Figure[l6]with each light curve being overlapped on the y magnitude. 
The spe ctrum became that for fre e-free emission about 4-5 days after the outburst while it was that for blackbody during the first 
3 days (Gallag her & Nev II1976I) . So, we shift each observed light curve down by /S.m\ and overlap it on the observed y light 
curve between day 6 and 30. Then, we obtain the color of Am\ = m y -m\ as listed in Figure PT6l 
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FIG. 20. — The V — R, V —I, V —J, V-H, and V — K (open triangles) and the y—R, y—I,y—J,y — H, and y — K (filled triangles) for the four classical novae: 
(a) VI 500 Cyg, (b) VI 668 Cyg, (c) VI 974 Cyg, and (d) GQ Mu s. Op en and filled circles represent, respectively, the intrinsic and the reddening corrected colors 
as expected from free-free emission (see Table|4). See equation (AlJ for the relation between free-free color and reddening). Note that the H and S'-band is the 
least contaminated by emission lines, so that the y—H and y—K colors are in a reasonable agreement with the theoretical values. See text for more details. 



TABLE 4 
Colors of free-free light curves 



color" 


intrinsic 


coefficient 


(m v -m x ) 


((M v -M x )o) 


(ca) 


V-V 


. 





V-R 


. 0.18 


0.7812 


V-I 


. 0.39 


1.6058 


V-J 


. 0.90 


2.2258 


V-H 


. 1.33 


2.5575 


V-K 


. 1.85 


2.7528 



' color is calculated from (my -m\) = {My -M\)o + 
C\E(B — V), where (My —M\)o is the intrinsic color of 
the free-free spectrum and c\ is the reddening coeffi- 
cient 

The intrinsic color indices, (My - M for V-R, V-I, V-J, V-H, and V — K of free-free emission (Table|4| are plotted in 
Figure l20h (open circles) together with E(B-V) = 0.45 (filled circles) calculated from equation (lAlb . 

The V-I, V-J, V-H, and V — K color indices obtained from Figure [16] are plotted in Figure l20k (o pen triangles). If strong 
emission lines are present in the spectrum, observed colors deviate from those calculated from equation (IA1 b - In fact, the V band 
is not emission-line free, being contaminated by strong emission lines even between day 6 and 30. The amount of contamination 
can be estimated from y— V = 0.3 because the y band is almost emission-line free. This value tells us that strong emission lines 
contribute about 30% of the energy flux to the V band. The y-I, y-J, y-H, and y—K color indices are also plotted in Figure 
|20t (filled triangles). 

A somewhat larger contamination is present in the / and J bands while the H and K bands are not s o hea vily contaminated by 
emission lines. For V1500 Cyg, the y—H and y—K colors are in a reasonable agreement with equation (IA1I) with E(B-V) = 0.45. 

V1668 Cyg 

Figu re [TTl sho ws y, V, I, J, H, and K li ght curves of VI 668 Cyg. Here the y magnitude obs ervations are taken from 
[Gallagher et al. ( 1980), the V data are from[Mallama & Skillman ( 1971) and iDuerbeck et al. I dl980l) . the / magnitudes from 
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iDeroux I ([1978), and the J, H, and K magnitudes from iGehrz et al. I d!980l) . Each light curve is shifted down to overlap with 
the y light curve between day 6 and 30 because the spectrum was no longer that for free-free emission after day 35 due to dust 
formation. The colors are summarized in this figure and are also plotted in Figure l20b . The / band is strongly, the J band is 
slightly, while the H and K bands are not heavily contaminated by emission lines. For VI 668 Cyg , the y — H and y—K color are 
in reasonable agreement with equation (lAlb together with E(B — V) = 0.40 (Stickl and et al. Ifl981h . Note also that the V band is 
contaminated by emission lines as suggested from y- V = 0.4. 

VI 974 Cyg 

For V 1 974 Cyg, no y and / photometr y is available. Figur e l20b shows the V — J, V—H, and V — K co lors obtained from Figur e 
[T8l where V magnitudes are taken from lChochol et alTI (1 19931) and the J, H, and K magnitudes are from lWoodward et al. Id 1997b. 
Open triangles denote the observed colors, which follow quite well the calculated ones with E(B — V) = 0.32 (Chochol et al. 
1993). However, the V band is contaminated by emission lines even at early phases as alre ady shown in VI 500 Cyg and V1668 
Cyg. If we subtract this contribution of AV ~ 0.3 (estimated from the spectrum in Fig. 2 of lBarger et al. I[l9 93), i.e., y—V = 0.3, 
we obtain the corrected colors (filled triangles) as shown in Figure l20b. The J band is h eavily contaminated by O I, Paschen /3, 
and 7, while the H and K bands are not so heavily contaminated (Woodw ard et al. Il 1997b . This explains why the y—H and y—K 
colors, corrected by AV =y—V= 0.3, are in good agreement with the expected values from free-free emission. 

GQMus 

In Figurel20h we plot (as open triangles) the observed V-/, V—J, V-H, and V—K colors from Figures[l9] and compare them 
with the intrinsic colors from free-free emission calculated from equation dAU for E(B-V) = 0.55 (filled circles). The observed 
colors are in reasonably good agreement with the calculated ones. 

However, the V ban d is not emission-line f ree, being contaminated by strong emission lines even in the early phase (e.g., 
~ 30% from Fig. 2 of Whit elock et al. Ifl984b . This contributes AV = y—V ~ 0.3 to the colors, so we subtract this from the 
original value of V, and obtain the corrected colors of y — /, y—J, y—H, and y—K (filled triangles) as shown in Figure l20b. The 
/ band is strongly co ntaminated by O I line and the J band is also very strongly contributed by Paschen j3 line (e.g., Fig. 3 of 
IWhitelock et al. 111984b . Consequently, both the y—I and y—J colors (corrected by AV =y — V= 0.3) are 0.8 and 1.3 mag below 
the calculated ones, respectively. On t he other hand, the H and K bands are not heavily contaminated by strong emission lines 
(e.g., Fig. 3 of lWhitelock et aini!984h . This explains why both the y—H and y—K colors corrected by AV = y— V = 0.3 are in 
good agreement with equation dAlb of E(B-V) = 0.55. 
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